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Abstract: As a simplest extension to the mass-term curvaton model, the curvaton model 
with a polynomial potential can relax the restricted constraint from PLANCK due to the 
non-linear dynamics of curvaton field before it decays. We find that there is still a big 
room for producing a large negative g^L, but not positive <7nl- For example, we only need 
around 10% "tuning" for — <7nl > 10 4 . 
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1. Introduction 

Even though inflation [|l| is an elegant model to explain the well-known puzzles in the 
hot big bang model, the mechanism for generating the observable anisotropies in cosmic 
microwave background radiation (CMBR) and forming the large-scale structure has not 
been well-established. In general, one may expect that there should be many light scalar 
fields (compared to the Hubble scale) during inflation and the quantum fluctuations of 
some of them finally seed the structure formation in our universe. Since the amplitude 
of quantum fluctuations of the light scalar field is proportional to the Hubble parameter 
which almost does not evolve during inflation, a nearly scale-invariant primordial density 
perturbation, or equivalently the spectral index of primordial curvature perturbation n s ~ 
1, can be taken as a strong prediction of inflation. This prediction has been confirmed by 
the CMB observations. For example, the full analysis [j^] of pre-PLANCK data, including 
the 9-year Wilkinson Microwave Anisotropy Probe (WMAP) |J, South Pole Telescope 
(SPT) 0], Atacama Cosmology Telescope (ACT) ||, Baryon Acoustic Oscillation ||] and 
Hq prior from HST project pj, implies 

n s = 0.961 ± 0.007 at 68% CL, (1.1) 

which is almost the same as that from PLANCK || || 

n s = 0.9603 ± 0.0073 at 68% CL. (1.2) 

Inflation also predicts gravitational wave perturbation whose amplitude compared to the 
scalar perturbations is measured by the so-called tensor-to-scalar ratio r. Unfortunately 
the gravitational wave perturbations have not been detected and PLANCK |j| sets an 
upper bound on it as follows 

r < 0.11 at 95% CL. (1.3) 

As the simplest setup of inflation, canonical single-field slow-roll inflation is governed 
by a canonical scalar field <f> (inflaton) and the expansion rate, the Hubble parameter, is 
determined by the potential energy of <p. During inflation inflaton slowly rolls down its flat 
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potential in order to achieve enough e-folding number for solving the puzzles in the hot big 
bang model. If the cosmic structure is completely originated by the quantum fluctuations 
of 0, it predicts that the curvature perturbation at the non-linear orders must be very 
small because large self-interaction of inflaton field implies a steep potential and breaks 
the slow-roll conditions ||10|| . Quantitatively a well-understood ansatz of non-Gaussianity 
has a local shape and the curvature perturbation C( x ) can be expanded to the non-linear 
orders as follows 

C(x) = Cl(x) + ~/nlCI(x) + ^nlCI(x) + • ■ ■ , (1.4) 

where Cl( x ) denotes the linear, Gaussian part of curvature perturbation. The sizes of the 
non-linear order perturbations are measured by the non-Gaussianity parameters, such as 
/nLj 5nl and so on. Different from single- field inflation, multi-held inflation model can 
easily produce large local non-Gaussianity. 

In this paper we focus on a well-known multi-field inflation model, namely curvaton 
model [11, [l3|, 14, 15] in which the final adiabatic curvature perturbation is generated 



by curvaton field a in the radiation dominant era far after the end of inflation. Usually a 
large local non-Gaussianity is expected in the curvaton model. Recently PLANCK data 
fl~6f| provides a stringent constraint on the local non-Gaussianity: 

/ NL = 2.7 ± 5.8 at 68% CL. (1.5) 

It implies that a large local bispectrum is unlikely. Even though the curvaton model can 
still fit PLANCK data well, there is no doubt that it has been tightly constrained. On the 
other hand, the size of trispectrum has not been constrained significantly. There are two 
shapes of local trispectrum which are measured by #nl and tnl respectively. Here we only 
consider the case in which the curvature perturbation is originated by single source and 
thus tnl = (|/nl) 2 - The constraint on tnl from PLANCK is 

r NL < 2800 at 68% CL (1.6) 

which is much looser compared to that from /nl- The constraints on ^nl is hopefully to 
be done in the near future. 

In S the curvaton model with quadratic potential was discussed. Here we consider a 



curvaton model with a polynomial potential which has been widely discussed in [14, 15, 17 



18, fil|, pO, 21]. In Sec. 2, we will investigate the constraint on such a model and then figure 



out the prediction of $nl in the constrained curvaton model. More discussion is given in 
Sec. 3. 

2. Qnl in the constrained curvaton model 

In this section we focus on the curvaton model with a polynomial potential, as the simplest 
extension to the model with quadratic potential, 

V(a) = I m V + Am 4 , (2.1) 



2 \m 
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here m is the mass of curvaton and A is the dimensionless coupling constant. Here m 2 is 
positive, but A can be positive or negative. For simplicity, the size of the self-interaction 
term compared to the mass term is characterized by 
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In this paper the subscript * denotes that the quantity is evaluated at the time of horizon 
exit of relevant perturbation mode during inflation. The equation of motion of curvaton 
field a after inflation is given by 
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The correction from the self- interaction term in the above equation is small if \s\ <C 2/n. 

Once the self-interaction term is taken into account, the non-Gaussianity parameters 
are significantly modified due to the non-linear evolution of curvaton field. In order to make 
our paper complete, we directly quote the results from [^Sf| . The amplitude of primordial 
scalar power spectrum generated by curvaton is 



r, 2 



and the non-Gaussianity parameters are given by 

/nl 

and 
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^l a ,D is the fraction of the curvaton energy density in the energy budget at the time of 
curvaton decay, 
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w(xq) + n(n — l)g(n, xq)s/2 
w(x ) + ng(n,x )s/2 
w(xp) + ng(n,x )s/2 
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(2.8) 

n{n-l){n-2)g{n,x )s/2, (2.9) 
n(n - l)(n - 2)(n - 3)5(71, x )s/2, (2.10) 



w(x ) = 2 1 / 4 r(5/4):r 1/4 J 1/4 (x ), 
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g(n,x )= ir2^' A T{b/4) n - l x 

x 
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and xo = wio = 1 denotes the time when curvaton starts to oscillate. Since < ^l a ,D < lj 
< f£> < 1. The spectral index of power spectrum generated by curvaton is 
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In the limit of s — > 0, the above results reduce to the model with quadratic potential 
and then 
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Considering m £ [0, 1], /nl > —5/4 and ^nl < 25/12. From the above two equations, we 
conclude that 
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In the curvaton model with quadratic potential c/nl is too small to be detected. 

From now on, we switch to the the case of s ^ 0. Usually we may expect that /nl 
should be quite large if td "C 1. However, from Eq. (fO|), /nl can be tuned to zero even 
for ro <C 1. In the limit of t\d — > 0, the corresponding value of s for /nl = shows up 
on the left panel of Fig. |]. Now one can check that hs + 3/i2 7^ which implies that <?nl 
can be arbitrarily large. The corresponding value of $nl x r £) is shown on the right panel 
of Fig. ||. It indicates that <?nl can be negative with large absolute value. However the 
fine-tuning for s is needed inevitably for obtaining a large $nl- 

Now let's go away from the fine-tuning point of s in Fig. |]. Since /nl has been tightly 
constrained by PLANCK data in Eq. ( |1.5| ), the value of m is constrained for different 
value of s. See the left panel of Fig. |2[ Even though /nl has been tightly constrained, 
there is still a big parameter space for the curvaton model with a polynomial potential. 
The allowed regions for ^nl are shown on the right panel of Fig. ^. From the right panel 
of Fig. ||, we find that the allowed value of $nl goes to infinity when s approaches to the 
fine-tuning point in Fig. |]. 

It is also interesting for us to investigate how much fine-tuning is needed for getting a 
large value of #nl- For example, the parameter space for — <7nl > 10 4 shows up in Fig. |3[ 
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Figure 1: The value of s and <7nl * rfj for /nl = in the limit of rjj — > 0. 

It implies that a large — #nl can be obtained by around 10% tuning around the fine-tuning 
point in Fig. |l]. 

Previously we mainly focus on the curvaton model with s<l which implies that the 
curvaton energy density is dominated by its mass term. One may expect that a large (?nl 
can be achieved if s > 1. However in such a case, /nl becomes large as well, and we have 
5nl ~ C(/nl)- See [|l5|, |2^, ^j] where the analytical and/or numerical calculations are 
presented in detail. Since /nl has been tightly constrained by Planck, we conclude that 
<?nl cannot be quite large in the curvaton model with dominant self-interaction term. 

Before closing this section, we also want to pay attention to the spectral index of power 
spectrum. From Fig. ||, the parameter s needs to be positive (around the fine-tuning point) 
in order to achieve a large value of ^nl, or numerically r] a /(m 2 /H%) ~ 0.8, 0.7, 0.6 for 
n = 4, 6, 8 respectively. It implies r) a > which makes the power spectrum bluer. Now 
one may worry about the constraint on the spectral index in Eqs. Ql . 1| ) and ( |1.2| ): a red 
tilted power spectrum is preferred at more than 5a level. In order to explain such a red 
tilted power spectrum, one need a large value of e which can be realized in the inflation 
model with potential U((p) ~ (jf . 1 In such an inflation model, e = where N is the 
number of e-folds before the end of inflation. For N = 50, e = 0.005p. For p = 4 and 
m/H* <C 0.1, the spectral index in the curvaton model is n s ~ 0.96 which can fit the data 
very well. 

3. Discussion 

In this paper we focus on the curvaton model with a polynomial potential. We find that the 
value of /nl can be tuned to zero as long as the curvaton self-interaction term has suitable 

1 In |i)J, the inflation model with p > 2 is disfavored at more than 95% CL because the predicted 
tensor-to-scalar ratio is much bigger than the bound in Eq. ( ]l.3| ). However we have to point out that the 
power spectrum is assumed to be completely generated by the quantum fluctuation of inflaton field <j> in 

|^| . But here the story is totally different: the power spectrum is assumed to be generated by the curvaton 
field. Denoting /3 = P C ,<r/P(, obs , r = 16(1 - /3)e < 0.1, even for p = 4, if (1 - #) < O(0. 1) which means 
that the curvature perturbation is mainly generated by curvaton field. On the other hand, this example 
tells us that some inflation models tightly constrained by PLANCK data might be relaxed in the curvaton 
scenario. 
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Figure 2: The allowed parameter space and prediction of ^nl in the constrained curvaton model. 



size compared to its mass term even when ro<l, and then ^nl can be arbitrarily large. 
However, in the curvaton model with dominant self-interaction term, this phenomenology 
does not happen and then <?nl cannot be quite large compared to /^ L . A fine-tuning to 
the strength of curvaton self-interaction is needed for obtaining a small /nl when < 1. 
Our numerical analysis implies that once such a fine-tuning is abandoned, <7nl can still 
be large if the strength of curvaton self-interaction is not far away from the fine-tuning 
point. We also notice that the parameter s must be positive at the fine-tuning point. It 
implies that the axion-type curvaton model cannot achieve a large <7nl when the constraint 
on /nl from Planck is considered. To summarize, even though the curvaton model with 
a polynomial potential has been constrained by PLANCK data, a large positive value of 
— <7nl can be generated without fine-tuning. 
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Figure 3: The allowed parameter space for — <jnl > 10 4 in the constrained curvaton model. 
Actually in many other models, such as the non-Gaussianity generated at the end of 



multi- field inflation model [22] and the general single- field ultra-slow-roll inflation model 
1 23] and some others in [24], ^nl is an independent parameter which is not related the 
value of /nl at all. In such kind of model, a large ^nl is still expected even though /nl 
has been tightly constrained to be around zero. 

Finally we want to point out that the canonical single-field slow-roll inflation predicts 
/nl = t|(1 — n s ) fpTofl . For n s = 0.9603, /nl = 0.0165 which is still far from the current 
sensitivity of detector. Checking this consistency is an important task in the future. 
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